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Abstract 

Contributions to the decay rates of P-wave charmonium states that are pro¬ 
portional to n,fa.g, where rif is the number of flavors of light quarks, are 
calculated in the framework of nonrelativistic QCD using the threshold ex¬ 
pansion method. Dimensional regularization is used to regularize the infrared 
divergences that arise from the emission of a soft gluon. Our results are 
consistent with the original calculations of Barbieri et al. 
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I. INTRODUCTION 


The nonrelativistic QCD (NRQCD) factorization formalism JIJ provides a systematic 
framework for analyzing annihilation decay rates and inclusive production rates of heavy 
quarkonium. These rates are factored into short-distance coefficients and long-distance 
NRQCD matrix elements. The short-distance coefficients can be calculated using perturba¬ 
tive QCD as a power series in the strong coupling constant a s (m c ) at the scale of the heavy 
quark mass. The matrix elements scale in a definite way with v, the typical relative velocity 
of the heavy quark. The NRQCD factorization formalism therefore organizes the decay rate 
or production rate into a double expansion in powers of a s and v. 

The threshold expansion method [U is a general method for calculating the short-distance 
coefficients which fully exploits the NRQCD factorization formalism. In this method, a 
quantity that is closely related to the creation or annihilation rate of a cc pair near threshold 
is calculated using perturbation theory in full QCD and then expanded in powers of the 
relative 3-momentum q of the c and c. Matrix elements of 4-fermion operators in NRQCD 
are also calculated using perturbation theory and expanded around the threshold q = 0. 
The short-distance coefficients in the factorization formula are then determined by matching 
these expansions in q order by order in a s . Finally, the NRQCD matrix elements for specific 
quarkonium states are simplified by using rotational symmetry, heavy-quark spin symmetry, 
and the vacuum-saturation approximation. 

In calculating the short-distance coefficients beyond leading order in a s , ultraviolet diver¬ 
gences and infrared divergences inevitably arise and need to be regularized. In perturbative 
calculations, the most convenient method for regularizing both ultraviolet and infrared di¬ 
vergences is dimensional regularization. We recently generalized the threshold expansion 
method to N spatial dimensions, so that dimensional regularization can be used consis¬ 
tently in quarkonium calculations [j3|]. The QCD side of the matching condition and the 
matrix elements on the NRQCD side are calculated in N dimensions using perturbation 
theory and then matched to obtain the short-distance coefficients. After renormalization of 
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coupling constants in QCD and NRQCD, the short-distance coefficients may have poles in 
N — 3, which must be removed by renormalization of the 4-fermion operators in NRQCD. 
One must take care to avoid simplifying the matrix elements of these operators using iden¬ 
tities that are specific to 3 dimensions until after these renormalizations have been carried 
out. We illustrated our method by calculating gluon fragmentation functions to S-wave and 
P-wave charmonium states, resolving a discrepancy between two previous calculations of the 
fragmentation function for g —> XcJ- 

There are also discrepancies in the literature between various calculations of the anni¬ 
hilation decay rates for P-wave states. In the original calculations by Barbieri et al. 
either the binding energy of the cc pair or the energy resolution of a gluon jet was used as the 
infrared cutoff. In recent calculations by Huang and Chao || and by Petrelli 0, dimensional 
regularization was used for the infrared cutoff. The discrepancies between their results and 
those of Barbieri et al. appear in the terms proportional to nga^, where rif is the number 
of light quark flavors. This is the term that is sensitive to the infrared cutoff. In this paper, 
we use the threshold expansion method in conjunction with dimensional regularization to 
calculate the term proportional to Ufa We find that the discrepancy between the results 
of Barbieri et al. and the recent calculations using dimensional regularization is due simply 
to different definitions of a color-octet NRQCD matrix element. 

The outline of the rest of this paper is as follows. In Section II, we review the threshold 
expansion method in N dimensions as it applies to decay rates. In Section III, we calcu¬ 
late the terms proportional to nja 2 in the short-distance coefficients of color-octet matrix 
elements and the terms proportional to and nga^ in the short-distance coefficients of 
color-singlet matrix elements. In Section IV, we apply our results to the decay rates of 
spin-singlet S-wave and spin-triplet P-wave states. Finally, we compare our results for the 
P-wave states with the previous calculations in the literature. 
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II. THRESHOLD EXPANSION METHOD 


The threshold expansion method for calculating the short-distance coefficients in in¬ 
clusive production cross-sections was developed in Ref. 0 and generalized to N spatial 
dimensions in Ref. [0|. In this section, we review this method as it applies to annihilation 
decay rates, since many of the formulas differ a little from the production case. 

The annihilation decay rate for the charmonium state H can be written in the factorized 
form JIJ 


r (H) = 


2 M, 




h ™ rn d c mn 


( 1 ) 


where Mh is the mass of the state H and d mn is the mass dimension of the operator O mn . 
The matrix elements (H\O mn \H) are expectation values in the quarkonium state H of local 
4-fermion operators that have the structure 


Omn = ^K-'lnX , 


( 2 ) 


where i/j and x are the field operators for the heavy quark and antiquark in NRQCD, and 
/C n and /C'l z are products of a color matrix (1 or T a ), a spin matrix, and a polynomial in the 
gauge covariant derivative D in N dimensions. The spin matrix is either the unit matrix 
or a polynomial in the Pauli matrices a 1 . The Pauli matrices in N dimensions satisfy the 
anticommutation relations {<jj, <jj} = 28 lv — In 3 dimensions, they also satisfy 

the commutation relations [a t , aj] = 2i i,j,k = 1,2,3. The commutation relations 

can be used together with the anticommutation relations to reduce all spin matrices to a 
linear combination of 1 and a 1 . However, since the commutation relations are specific to 
N = 3 dimensions, they should be used to simplify NRQCD matrix elements only after all 
poles in N — 3 have been removed from the short-distance coefficients. 

For the purpose of calculating short-distance coefficients, it is convenient to define the 
states | H) = |if(P = 0)) in (|l|) so that they have the standard relativistic normalization: 


H( P') H( P)) = 2 E P (2tt) n 8 N (P - P') , 


(3) 
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where Ep = y Mjj + P 2 . With this choice of normalization, the matrix elements in ([[]) dif¬ 
fer from the standard NRQCD matrix elements (H\Oi( 2S+1 Lj)\H) and (H\0 8 ( 2S+1 Lj)\H) 
introduced in Ref. jl|. The relation between them is discussed in Appendix B of [0. Up to 
corrections of relative order u 2 , the difference is a simple multiplicative factor. The renormal¬ 
ization of the operators O mn makes them depend on the renormalization scale /i of NRQCD, 
as indicated by the superscript (/i) on the matrix elements in ([!]). This dependence will be 
suppressed whenever it is not essential. 

The factors of m c in (i are chosen so that the short-distance coefficients C mn are di¬ 
mensionless. The dimension of the operator O mn in (||) is d mn = 2N + D, where D is the 
number of covariant derivatives in the operator. Since the coefficients C mn take into account 
the effects of short distances of order l/m c , they can be calculated as perturbation series in 
the QCD coupling constant a s (2m c ). The coefficients in the perturbation series depend on 
In (n/m c ) in such a way as to cancel the //-dependence of the matrix elements. 

The short-distance coefficients for annihilation decay rates can be determined by match¬ 
ing perturbative calculations of cc scattering amplitudes. Let cc(q, £, rf) represent a state 
that consists of a c and a c with spatial momenta ±q in the cc rest frame and spin and color 
states specified by the spinors £ and ■//. The standard relativistic normalization is 


c(qi,£) c (q2,h) 


4 E qi E q2 (2 tt) 2N 5 N (c[i - c( 1 )S N (q 2 ~ q' 2 ) , ( 4 ) 


where E q = y^m 2 + q 2 . The spinors are normalized so that = 1, and similarly for ■//, 
and r]'. Using the abbreviated notation cc = cc(q, £, 7/) and cc' = cc(q', £', r(\ the matching 
condition in the threshold expansion method of Ref. [|J is 


pQCD 


£(2? t) n+1 s n+ \p - k x ) {T c ^ x y 

X 

= E C d miP L (cc'W'K.'Lxx'lCnWcc)^ 

L — 'm LL mn iV - 1 


YYl LL mn 

mn "°c 


(5) 


pNRQCD ’ 

where P = (2E q ,0) is the momentum of the cc pair, T C c^x is the T-matrix element for 
annihilation of the cc into the state X consisting of light partons, and k x is the sum of the 
momenta of the outgoing partons. The sum over X on the left side of (|5|) includes integration 
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over the phase space of the final state partons and sum over their spin and color quantum 
numbers. By the optical theorem, the left side of (|5|) is proportional to the imaginary part 
of the annihilation contribution to the T-matrix element T^^cc'- Specifically, it is the sum 
of all cut diagrams for which the cut does not pass through any heavy quark lines. 

To carry out the matching procedure, the left side of (j|) is calculated using perturbation 
theory in full QCD, and then expanded in powers of q and q'. The matrix elements on the 
right side of (||) are calculated using perturbation theory in NRQCD, and then expanded in 
powers of q and q'. The short-distance coefficients C mn are obtained by matching the terms 
in the expansions in q and q' order by order in a s . 

The calculation of the short-distance coefficients can be simplified by averaging both 
sides of the matching condition (|j) over rotations of the vectors q' and q and the spinors £, 
7], and rf that specify the states of the initial and final cc pairs. On the NRQCD side of 
the matching equation, the average over rotations can be accomplished simply by restricting 
the operators O mn to be rotationally invariant. 

In the perturbative calculations of the matching condition (|5|), infrared and ultraviolet 
divergences can appear on both sides of the equation. Since NRQCD is constructed to be 
equivalent to full QCD at low momenta, the infrared divergences on both sides must match. 
They therefore cancel in the short-distance coefficients C mn . Any ultraviolet divergences on 
the left side are eliminated by renormalization of the QCD coupling constant and the heavy 
quark mass. On the right side, the ultraviolet divergences are eliminated by renormalization 
of the gauge coupling constant and other parameters in the NRQCD lagrangian and by 
renormalization of the 4-fermion operators of NRQCD. 

The matching calculations are particularly simple if dimensional regularization is used 
to regulate both infrared and ultraviolet divergences. Radiative corrections to the NRQCD 
matrix elements vanish identically, because after expanding the integrand of the radiative 
correction in powers of q and q', there is no momentum scale in the dimensionally regularized 
integral. The radiative corrections to the matrix elements do include infrared poles in 
e = (3 — N)/2 that match the infrared divergence on the QCD side of the matching condition, 
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but they are canceled by ultraviolet poles in e. Thus the only noncanceling contributions 
to the NRQCD side of the matching condition are the tree-level contributions of the matrix 
elements, including those matrix elements that arise from counterterms associated with 
operator renormalization. 

III. PERTURBATIVE MATCHING OF QCD AND NRQCD 

In this section, we use the threshold expansion method to calculate selected terms in the 
short-distance coefficients in the factorization formula We calculate the term propor¬ 
tional to nfd 2 s from the color-octet annihilation processes cc —> qq and the terms propor¬ 
tional to a 2 and rifa^ from the color-singlet annihilation processes cc —> qqg and cc —> gg. 
Dimensional regularization is used as a cutoff for both infrared and ultraviolet divergences. 


A. Color-octet terms from cc —> qq 


The terms on the QCD side of the matching condition ([5]) that are proportional to Ufa 2 
come from the annihilation process cc —> qq. For these terms, the left side of (^]) reduces to 


[ [ (2vr ) N+1 5 N+ \P - l^iTc-c'^YTc-^ , 

J i\ 2 


( 6 ) 


where l\ and li are the momenta of the q and q and the remaining sum is over their colors, 
spins, and flavors. The symbol f k denotes the N -dimensional integral over the phase space 
associated with the momentum k: 



(7) 


The term in (|S|) that is proportional to Ufa 2 is given by the cut Feynman diagram in 
Figure 1. The T-matrix element for c{p)c(p) —> g(Zi)g(7 2 ) is 
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The c and c have 4-momenta p = (E q , q) and p = (E q , — q) where E q = \Jm 2 c + q 2 . The 
coupling constant in (D) has been written g s p, where e = (3 — N)/2 and // is the scale of 
dimensional regularization, so that g s remains dimensionless in N dimensions. The integrals 
over £i and £ 2 in (b]) can be carried out by using the energy-momentum-conserving delta 
function: 


/ f (2 7 r) iV+1 h 7V+1 (P-/ 1 -/ 2 )^h(/ 2 )7 v r fe n(/ 1 )h(/ 1 )7 /i T a n(/ 2 ) 


= n f 


(N — 1) T(|) (P 


2 \ - £ 


167T 


(- P 2 g > + P,P t ) 5 ab , 


(9) 


8nNT(f) 

where P = (2 E q , 0 ). We have expressed this in Lorentz-invariant form for later use. Inserting 
this into the QCD side of the matching condition, it reduces to 

1 


(N — 1) T(|) 4 4e 
n f TPPrPk 9 s V 


/ E p ~ f 


4 ? 


E l 


u(p)YT a v(p)v(p)^ tl T a u(p) . 


( 10 ) 


32ttAT (f) 

Using the formulas in Appendix A, v(p)'y fl T a u(p) can be expressed in terms of nonrela- 
tivistic Pauli spinors. Expanding to linear order in q, we obtain 


v(p)^T a u{p) « —2 m c g' 11 p . 


( 11 ) 


The QCD side of the matching condition ( lOj) then reduces to 


2ir (JV — 1) r(|) 


2 / 2 4e 


n f - - ,TT,,N^ ~ a sP" l ) C'o-T a r]' ■ if crT a f . 


NT( 

Setting N = 3, this simplifies to 


'N\ 
2 > 


47T 


( 12 ) 


a l poVrt ■ i,'oT“£ . 


(13) 


Dividing by (2 m c ) 2 to account for the difference in the normalization of states, we reproduce 
twice the leading term in ImAd given in (A13) of Ref. JI[]. 


B. Color-singlet terms from cc — ► qqg 

Terms proportional to n/cxf on the QCD side of the matching condition come from the 
processes cc —> qqg and cc —> qq. For the process cc —> qqg, the left side of (|5|) can be 
written 
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(14) 


[ [ [ (2t T ) N+1 5 N+1 (P-k-h-£ 2 ) 
J k J J in 


cc *qqg > 


where £i, £ 2 , and k are the momenta of the q, q, and gluon, respectively, and the remaining 
sum is over the colors, spins, and flavors of these partons. The term proportional to nfa^ 
comes from the cut Feynman diagrams in Figure 2. The T-matrix element in Feynman 
gauge for c(p)c(p) —> q(h)q(h)g(k) can be written 


Tcc^qqg = A^ g * g ( " 2 (g 8 lfu(h)'y v F a v(l 2 ) 


(£i + £i 


(15) 


where A^_> g * g is the amplitude for the process c(p) c(p) —> g*(l)g(k): 
^-n 7*9 = 3 s V e £(k) v(p) 


~ P + m c )Y^ rphrpa 


2 p ■ k 


7 M (jf — ft + m c ) 7 " 


+ ' '' 2f~k T ° Tt “ (P) ' 


(16) 


We will compute the matching condition for the case of color-singlet cc pairs. The color 
matrices T a T b and T b T a in ( |T 6 | ) can then be replaced by 5 ab / 6 . We will further simplify the 
matching condition by averaging the initial state cc(q, £, 7 ) and the final state cc(q', ?/) 
over the rotation group. 

The integral over £\ and 1 2 in (|T|) can be carried out by inserting the identity 




(17) 


After using ([|) with P replaced by i to integrate over £\ and £ 2 , we obtain 
' 3 ^ ' dJL 2 




x <Fy(i^) a®) 

The remaining sum is over the colors and spins of the real gluon. The £ fl £ 1/ term gives no 
contribution due to current conservation. 

We proceed to carry out the nonrelativistic expansion of the factor in ([18]) that involves 
the amplitude A^_, g * g . Using the identities in Appendix A, we expand ( 0 ) to linear order 
in q: 
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A 

Slcc^g*g 


12m c ko 


e a u (k){ m c g^g" 3 rf{[a\a J ],k ■ cr}£ 


■j- k { (k»g vj + g^k u - 2m cfi M y° - k j g' 11 ') 

- k\g^g^ - <fV) + F(<f 0 <T - g^g*) 


k 0 g^g uj + (2 m c -k 0 )g”g' 


rfq'o 3 ^ 


( 19 ) 


In N > 3 dimensions, the spin matrix {[<7®, cr 7 ], cr k }, which is totally antisymmetric in its 
three indices, is linearly independent of 1 and cr\ i = 1,..., N. In 3 dimensions, it reduces 
to the unit matrix multiplied by —4 ie l ^ k . 

After multiplying A^^ g * g by (A^^ g * g )*, the spinor factors can be simplified by averaging 
the vectors q and q' and the spinors £, g, and g' over the rotation group. The factor 
£'V fc a l g'g'q l a^ reduces to a linear combination of three independent spinor factors: 

1 


£/f q'k(j l rfrp q l (J^C = 


N(N — 1)(N + 2) 

x ((AT + 1 )6 ij 5 kl - 8 ik 8 jl - 5 il 5 jk ) £' f q' • ag'rfq ■ cr£ 
+ ((IV + 1 )8 ik 8 jl - 8 ij 8 kl - S il 6 jk ) ^q ,m a n ri'rfq m a n t 
+ ((A" + 1 )8 il 8 jk - 8 ij 8 kl - 8 ik 8 jl ) i'^ m a n g'g ] q n a m i 


( 20 ) 


The spinor factors £' ji q' l a m rj'rf{[(j\a : ’],a k }£, and ^{[a l ,a m ],a n }g'g^q l cr^ average to zero, 
while {[cr l , cr m ], a n }r]'rf{[cr\ pi], cr k }t; reduces to a single spinor factor: 


£4{[cb, cr m ] 5 a n }g'g^ {[a 1 , cr- 7 ], cr k }t; 

8 il 8 jl 


N(N — 1)(N — 2) 


8 kl 

^im fijm km 

8 in <p n 8 kn 




( 21 ) 


In 3 space dimensions, the spinor factor on the right side of ([21]) can be simplified by using 
the commutation relations of the Pauli matrices: 


£'t{K,«T = -geeyye n = 3 . 


( 22 ) 
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However this simplification can only be made if there are no poles in N — 3 multiplying the 
spinor factor. After averaging over the rotation group using (]20j) and (0) and simplifying 
the Lorentz algebra, we obtain 




E (AS. 


tfia 

>9*9 J >A *' CC— >9*3 


1287T 2 2 4e ^ 

9N(N + 2) asfi x 2 


x ( 6(n + 2)x 2 


4 (IV + 2)x + (N — l)x 2 —— ■ CTTj'rjhi ■ cr £ 


m. 


±{N + 2)(N- 1)(1 - x) - (2N 2 + N- 9)x 2 


mi 


^q /m a n V Wq m a n Z 


4(TV + 2)(N- 2)x - (2N 2 -N- 7)x 2 ] — ^q ,m a n r) , rfq n a m £ , (23) 


mi 


where x = k 0 /m c . 

At this point, the integrand of (|18|) has been reduced to a function of £ 2 = 4 E q (E q — k 0 ). 
We can therefore carry out the phase space integral over k and £: 




N—3 


p2 


(24) 


Changing the integration variable to x = 1 — £ 2 /(4 E 2 ) and taking the nonrelativistic limit, 
the QCD side of the matching condition reduces to 

(n- i)r 2 (- 3N ' 9X_2e 1 


n f 


32n 2 NT 2 (f 


1L a u 2e 

jV\ 


m r 


\—1—e 


dxx 1 2e (1 — x)~ 

x(-v)EbS-,-,)‘AT 


■9 9 


(25) 


Infrared divergences arise from the endpoints of the integral in (^5|). The infrared di¬ 
vergences from x —> 1 arise when the q and q are collincar with the virtual gluon, which 
is almost on-shell. They are canceled by collinear infrared divergences from radiative cor¬ 
rection to the process cc —> gg, which will be calculated in subsection IIIC. The infrared 
divergence from x —> 0, which appears only in the coefficient of ^q lrn a n 'g l rfq rn a n £ t) arises 
when the emitted real gluon is soft. It must be matched by an infrared divergence on the 
NRQCD side of the matching condition. The analytic expressions for the integrals over x is 
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dx x 


n— 1—2e 


(1 -x) 


—1—e 


r(w - 2 e)r(—e) 

T(n — 3e) 


(26) 


Integrating over x and keeping only the terms that survive in the limit e —> 0, we obtain 


n f n,r>^T , oWv, 


4(iv-i)r 2 (|) , ^ t m: 

9N 2 (N + 2)r 2 (^) 


4(/r 


x- 


1 

tin 


- 6(5 + 3e) cr a ], (J t }g'rf{[cj r , ^1, ^K) 


- 2(11 - 4e) — £ ,f q' • ar)'rfq ■ <r£ 


nr 


2(16 - 21 e)—^q' m a n r)'rfq m a n £ 


m- 


2(6 - 17e) 


nr 


(27) 


The subscript IR on the pole in e is a reminder that it has an infrared origin. 


C. Color-singlet terms from cc —> gg 


The remaining terms on the QCD side of the matching condition that are proportional 
to rifoi i s come from the process cc —> gg. For this process, the QCD side of the matching 
condition (|5]) reduces to 

5 jj{ 2v) N+1 f> N+ \ p - k - . ( 28 ) 

where k and £ are the momenta of the gluons, the sum is over their spins and colors, and 
the factor of | comes from Bose statistics. 

We first consider the QCD side of the matching condition at leading order in a s . At 
this order, it is given by the cut Feynman diagrams in Figure 3. The T-matrix element for 
c(p)c(p) -> g(k)g(l) is 


T- = A^ a c a (£) 

1 cc-^gg s\ C c—>gg J 


(29) 


where Acz^ gg is identical to (0) except that £ 2 = 0. The nonrelativistic expansion of the 
matching condition is similar to but simpler than that for cc —> qqg, which was carried out 
in subsection IIIB. The final result is 
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8*r(i) 


. 2. ,4e 


m: 


9iV(Ar + 2)r(f) v 47r / 

r 6 (JV + 2) (- { [V, V], a ‘} VV { K, V], <t‘ K) 


+ (5JV + 7) —j f' q' ■ (TijVq ■ V 
m z c 

+ (2AT 2 + AT - 9) _L ^t ? ^ (T n r7 / ?7 t( g m (7 n + 

mi 


( 30 ) 


Setting N = 3, this simplifies to 
167T 


135 


a c 


15 + 11 —- £ 4 q' • trrfr^ q ■ <x£ 




+ 6 ^q' m a n r ] l r 1 \q m a n + q n a m )Z 


mi 


(31) 


Dividing by (2 m c ) 2 to account for the difference in the normalization of states, we reproduce 
twice the corresponding terms in ImAd given by the sum of (A9a) and (A9b) of Ref. [jjlj. 

The terms in (^) proportional to rifOp s come from the renormalization of the coupling 
constant in (J30J and from the cuts of the Feynman diagrams in Figure 2 that pass through 
2 gluon lines. The renormalization contribution in the MS prescription is obtained by 
multiplying (|30|) by Z*, where 


Z g = 1 + ^-a s (J— + hi(47r) - 7) ■ (32) 

The subscript UV on the pole in e is a reminder that it has an ultraviolet origin. In the 
diagrams obtained by cutting two gluon lines in Figure 2, the quark loops have ultraviolet 
divergences and collinear infrared divergences. With dimensional regularization, there is 
a calcellation between the ultraviolet and infrared divergences and the diagrams vanish. 
Making this cancellation explicit, the sum of the diagrams can be expressed as ([fCj) multiplied 
by the factor 


(---)■ 

5vr \e uv €irJ 

The ultraviolet pole in e is cancelled by the renormalization factor (|5^). 
contributions proportional to n/cxf, we obtain 


(33) 


Adding the two 
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n f 


8r (l 


27iV(iV + 2)r( 


N\ a s^ 


3,4 e(a\Y (± 
4vr/ \e IR 


+ ln(47r) — 7^ 


x 


6(iv + 2 ) 


+ (5 N + 7) -7 f' q' ■ ar/'fq ■ <r£ 


mt 


+ (2iV 2 + jV - 9) — ^q ,m a n r]'rf(q m a n + <f a m )£ 


mt 


( 34 ) 


On the QCD side of the matching condition, the total contribution proportional 
to 7ifa 3 is the sum of (^7|) and (|34]). The pole in e in front of the spinor factor 
£ t ^{[(T r 1 cr s ],cr t }ri , rf{[cr r ,cr s ] 1 cr t }t; cancels in the sum. The spinor factor can therefore be 
simplified using (|22|) . Adding (f27l) and (0) and keeping only the terms that survive in the 
limit e —> 0, we obtain 


n f 


sr(f) 


27N(N + 2)r((4) 


N\ a s h 


3, ,4e ( m l\ 


47T 


60(lll £: + IA'Wf 


■ 441 ln + 1) 5 ' ,q ' ■ CT,/ ’ ,tq ' ^ 


+ 


-20 


(± +H47l) - 1 y 6i (i n jL- + i 


mt 


^q' m a n ri'rfq m a n Z 


~ 24,>n £r c + 5) 


(35) 


D. Short-distance coefficients 


We now calculate the matrix elements that must appear on the NRQCD side of the 
matching condition (|5|) and extract their short-distance coefficients. We will find that the 

terms that are required for matching are 

ri&S) /'Aids') 

— (cc'I^XxVlcc) + jV _ 1 (cc'\^(rT a x ■ x^T a 7jj\cc) M 


m 


N- 


mt 


+ 


+ 


+ 


Ci 


(L 3 p) 


m 


N+l 


■(cc'|V’ t (-fD • a)xxK~0 ■ °)x\cc) 


2 „ + i xH-lD) m a n ip\cc) 

m c 

r (1,3p) «_> 

3 {ce’Wf-lDyWxx'HDra^cc) + 


m 


N+l 


(36) 
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We have placed a superscript (//) on the color-octet matrix element in anticipation of the 
fact that it will acquire dependence on the NRQCD renormalization scale through renor¬ 
malization. 


We consider the matching condition first for color-octet cc pairs and then for color-singlet 
cc pairs. In the case of color-octet cc pairs, the contribution proportional to nfO? s on the 
QCD side of the matching condition is given in (0). The spinor factor in (|l^) can be 
identified as the expression at leading order in a s and expanded to linear order in q and q' 
of the following NRQCD matrix element: 


(cc'\^crT a x ■ / \'^crT a '0|cc) ~ 4 m 2 c ^uT a r /' ■ rfcrT a £ t . 


(37) 


The tree-level expression for the matrix element in ([37]) is represented diagramatically in 
Fig. 4, with the dot representing the operator. Comparing (|T2D and (|36|), we can read off 
the short-distance coefficient of the matrix element in (]37|): 

' 4:7T /i 4 \ C 


= (: 
1 2AT(f) s V 


m„ 


(38) 


Taking the limit N 3, this reduces to 




(39) 


Since the short-distance coefficient is sensitive only to momenta of order m c , we have set 
the scale of the running coupling constant to /i = 2 m c . 

In the case of color-singlet cc pairs, the contributions proportional to a% and n/a^ on 
the QCD side of the matching condition are given in (|3l|) and (p5|) , respectively. The spinor 
factors in these expressions can be identified as tree-level expressions for NRQCD matrix 
elements, expanded to linear order in q and q'. The S-wave spinor factor is proportional to 


(cc'l^'x xVl cc) ~ 4m^ ^ rfr) 


(40) 


Comparing the sum of ( |5T| ) and (^3|) with (]3^), we can read off the short-distance coefficient 
of the matrix element in 


15 







c a ' s) = 


1 + 


- — In — 

3 2 m r 


CXs 

9j nf V 


(41) 


Choosing the QCD renormalization scale to be fi = 2 m c , this reduces to 

47T 0 ,_ s C 8n/ a s 


C°-' s) = f a’(2mj (l 


9 7T 


(42) 


The P-wave spinor factors in ([31]) and (|35|) can be identified with the matrix elements 


(cc' |'0 t (-|D • cr)x X f (-fD • cr)^|cc) ss 4m^ £ /f q' • crr/i^q • cr£, 
(ce'|^ t (-|B) m a n x X t (“|S) m cr>|cc) « Ami &qt m <T n rfrfq m <T n Z, 
{ cc'\^(-0) m a n X x f (“|B) n c7 m ^|cc) « 4m^ ^q ,m a n r]'^q n a m i . 


(43a) 

(43b) 

(43c) 


For the matrix elements ( |43a| ) and (43c), we can immediately read off the short-distance 
coefficients by comparing the sum of ([31]) and ([35]) with (]36|) . Setting N — 3 and /i = 2 m c , 
we obtain 


c-yv, = 44x a 2 (2mc) 0 . 58 n/ a 


135 

87T 


C , 3 i ’ 3P) = ^«s( 2m c) (l - 


99 7T 
3 7T 


(44) 

(45) 


The coefficient of the matrix element in ( |43a|) is more complicated because the coefficient 
of £' f q ' 7n o ri ri'rfq m a n ^ in (|35|) contains an infrared pole in e, indicating that it is sensitive to 
long-distance effects. Since an infrared divergence cannot appear in a short-distance coeffi¬ 
cient, that divergence must be matched by an infrared divergence in a radiative correction 
to a matrix element on the NRQCD side of the matching condition. Since the divergence in 
(|35|) has a coefficient proportional to the infrared-divergent NRQCD matrix element 

must have a short-distance coefficient proportional to rifOt 2 s . The only such matrix element 
is (cc'|^<rT a x • x^<TT a ^|cc), whose short-distance coefficient has already been determined 
in (|38|). Thus the infrared divergence on the NRQCD side of the matching condition must 
come from that term. 

If the initial and final cc pairs are in color-singlet states, the tree level expression © 
for the matrix element vanishes and the leading contribution comes instead from radiative 
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corrections. The leading contributions to the matrix element are represented by the diagrams 
in Figure 5. The expression for diagram 5a is 


V/i 2e ^a n T a T b r}'rfa n T b T a £ 


d N k 


(2n) N 2k 


• q' - q • k k ■ q' 


x 


E q — k — (q — k) 2 / (2m c ) + ie E q — k — ( q' — k) 2 / (2 m c ) + ie 


(46) 


where E q = q 2 /(2 m c ) = (q') 2 /(2 , m c ). Since the initial and final cc pairs are both color- 
singlets, we can replace T a T b and T b T a by 5 ab /6. The proper way to evaluate the diagram 
is to regard |q|, |q'| , and |k| to all be much smaller than m c . We must therefore expand 
out the denominators of ([K]) in powers of q /m Cl c(/m c , and k jm c before integrating over 
k. Keeping only terms up to linear order in q jm c and q '/m c , the diagram reduces to 


167T 
~ 9 ~ 


a s n 2e ^a n rjrfa n i 


d N k q • q' - q ■ k k • q' 


(47) 


(2ir) N k 3 

The integral is both ultraviolet and infrared divergent. It vanishes in dimensional regular¬ 
ization due to a cancellation between an ultraviolet pole in e and an infrared pole. Making 
these poles explicit, the diagram can be written 


a, (- --) (' 9”V‘? . (48) 

277r \euv e iR' 

The subscripts UV and IR on e indicate whether the pole is of ultraviolet or infrared 
origin. We have set N = 3 in the prefactor, since any finite terms obtained by expanding 
the prefactor in powers of e will cancel. The other 3 diagrams in Figure 5 give identical 
contributions. The final result for the matrix element is 

(cc'|V'Vr*x • xVrvicc) = ^(- -I)A'V,V?VC (49) 

2Y7T \€uv tiRS 

After multiplying by /m^ -1 , where is given in (|38|), we find that the infrared 

pole in e matches the one on the QCD side of the matching condition, which is given in (p5|). 

The ultraviolet pole in e in ([49|) indicates that the matrix element is ultraviolet divergent 
and therefore requires renormalization. It is the renormalized matrix element that appears 
in the NRQCD side of the matching condition (|55|). In the MS renormalization scheme, the 
relation between matrix elements of the bare operator and renormalized operators is 
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(cc'\^crT a x • ’x}crT a il)\ce) = n 4e ( (c^\^crT a x ■ x^crT^lcc)^ 


8 a, 


+ ln(4vr) - 7 ) (cc / |^ t (-|B) m a n xX t (-^) m a>|cc) 


O) 


(50) 


27Trml \e uv 

The superscripts (/i) on the matrix elements on the right side indicate that they are renor¬ 
malized matrix elements with renormalization scale /i. We will suppress this superscript on 
color-singlet matrix elements, since they do not require any renormalization at this order in 
a s . The fermion field operators in the bare matrix element on the left side of (|50|) have di¬ 
mension N/2. The fermion held operators in the renormalized matrix elements on the right 
side have dimension 3/2. The factor of n~ 4e on the right side of (^0j) compensates for the dif¬ 
ference between the dimensions of the two sides. Solving (|50|) for (cc']^crT a x • x^<rT a |cc) ^ 
and using (|49| ) and (|43b|) , we find that the renormalized matrix element is 


{cc'\^crT a x ■ x^<tT“^|cc)^ = 


32 ( 1 

-a s I — 


+ ln( 47 r) — 7 ^ ^q'" l a n rfif q m cr n £, ■ (51) 


277r \e IR 

Multiplying by C&V/m*- 1 , where C is given in fl38|) and keeping all terms that 
survive in the limit e —> 0, we hnd that the contribution to the NRQCD side of the matching 
condition is 


C& 3 s) 


m 


N— 


- {cc'\^crT a x • / \' t crT a '0|cc)^' ) = rif 


8T( 


x 


27N(N + 2)T(y 
1 


.2\ 

jvX/* 1 73 


m r 

47T 


1 7 

-20 (-K ln( 47 r) — 7 - 

^IR 5 


rnt 




(52) 

which 


We have expressed this in a form that makes it as easy as possible to match with 
is the term proportional to rifCe^ on the QCD side of the matching condition. In particular, 
we have expanded out a factor of (N — l)(iV + 2) to get the constant — | under the pole in 
e. We see that the infrared poles in e in (|35j) and (|52|) match. The remainder of that term 
must be matched by the term proportional to C/r’ P ' 1 in (|id|) . The resulting expression for 
the coefficient is 


riQ-’ 3p ) — 2/n 

^2 — ~^ a sW 


1 + 

9 2m c 27/ 7 r 


(53) 
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Shifting the renormalization scale of the QCD coupling constant to ji — 2 m c we obtain our 
final result: 


C , 2 i ’ 3P V) = 


The remaining logarithms of /i represent 
renormalization scale /i of NRQCD. 


10 j /i 49 \ ck/ 

9 2m c 27/ 7T 

the dependence of the matrix element on the 




IV. ANNIHILATION DECAY RATES 


The factorization formula (]T|) for the annihilation decay rate of a quarkonium state H, 
including all terms whose coefficients have been computed explicitly in Section III, is 
1 


T{H) = 


2 M h V ml 


(7 (Ids) /yd, 3 5) 

(H\^x xH\H) + - 3- {H\tpoT*x ■ x^ry\H)^ 


nit 


Cl 


(1, 3 P) 


nr 




+ 


c. 


a ?p) 




a 


m2 

(i, 3 p) 


(H\^(-0ra n x x\-0) m v n m) 


(H\^(-m m a n xx\- L 2 D) n a m m) + 


nr. 


(55) 


We have calculated the terms proportional to nfOi 2 s in C^-’ 3 ^ and the terms proportional 
to a 2 or nfOil in the other coefficients in ([55]). The ^-dependence of the coefficient C 2 ~’ P ' ) 
cancels that of the renormalized matrix element, which satisfies 


= xU _io ra n m} . (56) 

277T ml 

The relative importance of the various terms in the annihilation decay rates (|55|) depends 
on the quarkonium state. The magnitude of a particular term is determined by the order in 
a s of its short-distance coefficient and by the scaling of the matrix element with v, which is 
given by the velocity-scaling rules of NRQCD [1|]. Below, we apply this general formula to 
the annihilation decay rates of spin-singlet S-wave states and spin-triplet P-wave states. 
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A. Spin-singlet S-wave states 


The dominant Fock state of the g c consists of a cc pair in a color-singlet 1 S'o state. The 
largest matrix element that appears in the NRQCD factorization formula for the decay rate 
is therefore {gdi*'XX'rflVc) i which scales as v 3 . The next most important matrix element 
is {rjcK^XX 1 D 2 ^ + h.c.)|7/ c ), which is suppressed by a factor of v 2 . Matrix elements whose 
dominant contributions come from higher Fock states are suppressed by v 3 or more. For 
example, the matrix element whose dominant contribution comes from the ccg Fock state 
in which the cc pair is in a color-octet l P\ state is suppressed by v 4 , with one factor of v 2 
arising from the probability of the ccg Fock state and the other arising from the derivatives 
in the P-wave operator. The matrix element (r) c \^crT a x • / \'^crT a '0|?/ c ), whose dominant 
contribution comes from the ccg Fock state in which the cc pair is in a color-octet 3 S\ state, 
is suppressed by v 3 from the probability of that ccg Fock state. Thus, up to corrections of 
relative order v 2 , the annihilation decay rate can be written 

1 (7 (Ids’) 

T (Vc) = -— (Vc\^XX^\Vc) • (57) 

m c 

The short-distance coefficient, including the next-to-leading order correction proportional 
to Ufa 3 , is given in (^). The Ufa 3 term agrees with the complete next-to-leading order 
correction calculated by Barbieri et al. [£| and by Hagiwara et al. [0J. 

The difference between the matrix element in (0) and the standard NRQCD matrix 
element (r 7 c |(P 1 ( 1 S'o)| 7 c ) introduced in Ref. [|j] is discussed in Appendix B of fl2[. Up to 
corrections of relative order v 2 , the difference is simply an overall normalization factor: 

(Vcl^XxtylVc) ~ 4m c (rjclOiCSo)]^) ■ (58) 

Using the vacuum-saturation approximation, this matrix element can be related to the radial 
wavefunction of the g c evaluated at the origin: 

< I , c |CM 1 S 1 ,)k} « T|ft,.(0)| 2 . (59) 

This approximation is accurate up to corrections of relative order v 4 . 
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B. Spin-triplet P-wave states 


The dominant Fock state of the Xcj consists of a cc pair in a color-singlet i Pj state. This 
Fock state gives the dominant contributions to the color-singlet P-wave matrix elements in 
(|55|). Since the derivatives in the P-wave operators give a suppression by v 2 , the color-singlet 
P-wave matrix elements scale as v 5 . The matrix element a 'T a x ■ x^<rT a, 0|x c j), whose 

dominant contribution comes from the ccg Fock state in which the cc pair is in a color-octet 
3 S'i state, is suppressed by a factor of v 2 from the probability of that ccg Fock State. Thus, 
it also scales like u 5 , like the color-singlet P-wave matrix elements. All other matrix elements 
are suppressed by v 2 or more. Thus, up to corrections of relative order v 2 , the decay rate is 
given by the sum of the color-octet S-wave term and the three color-singlet P-wave terms in 


the NRQCD factorization formula (|55|) . 

We will use rotational symmetry and the approximate heavy-quark spin symmetry of 
NRQCD to show that the color-octet S-wave matrix elements (xcj\^ <J T a X ' X ] crT a ^\xcj) 
for the three XcJ states are equal up to corrections of relative order v 2 . Spin symmetry 
implies that the state Xcj(jz ), which is an eigenstate of J 2 and J z , can be expressed in the 
form 


Xcj{jz 


ls z\Jjz) Xc{lzS z )) , 


IzSz 


(60) 


where Xc{l z s z ) is an eigenstate of L z and S z . Thus the color-octet 3 Si matrix element can 
be written 


{Xcj\^<rT a x ■ X^T^lxcj) = 


2J + 1 


Y,^j{jz)\^ ] (TT a x ■ X f crT a ^\xcj(jz)) 


Jz 


2J + 1 


H H 'EiJjz l 1 ^; i s' z )(H z ; is z \Jj z ) 

jz l' z s' z IzSz 

x (xc(l' z s'z) v l T a XX ] ° l T a ^ XcilzSzj) ■ (61) 


Spin symmetry also implies that the matrix element on the right side of (^Tj) is proportional 
to Uj s , U Sz i, where U m i is the unitary 3x3 matrix that transforms vectors from the Cartesian 
basis to the spherical basis. Finally, rotational symmetry implies that matrix elements on 
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the right side of ([H|) must be proportional to Si/ z i z - Using the orthogonality relations of the 
Clebsch-Gordan coefficients, the equation (H) can be reduced to 

(XcjWcrT a x ■ X^Ty\xcj) ~ <Xco| ^<rT a X' X^T a ip\x c o) ■ (62) 


These relations hold up to corrections of relative order v 2 . 

We next show that the color-singlet P-wave matrix elements in (|55| ) can be reduced to 

• ^^ ^^ 

a single independent matrix element, which we choose to be (Xcol^' (—§D • (t)xX'(^ fD • 
cr )'01 Xco) - Using the expression (pD|) for the XcJ states, we can write 




2J+ 1 

* (xc(l z s 2 


EEE<4|u:;i4)(^;k-l4) 

jz l' z s' z l z s z 

xAs, 


(63) 


The vacuum-saturation approximation, which is accurate up to corrections of relative order 
v 4 can be used to express the matrix element on the right side of (|63]) as the product of 
{Xc{l , z s 'z)\x\~\D) m a n 'il>\Q) and (0\^(—^Dya j x\Xc(Ls z ))- Rotational symmetry and spin 
symmetry imply that these two matrix elements are proportional to Uj ls , and Ui z iU Sz j, 
respectively. Thus the tensorial structure of the matrix element ( |63| ) is completely deter¬ 
mined. The proportionality constant can be deduced by taking the special case J = 0, i = j 
and m = n, summed over i and m. The resulting formula is 


(XcM\-{0) i a>xx'(-kDra"iP\Xcj) ~ tt {xojW'HB' ' ‘BV’IXco) 


X 


2J + 1 


E E E«.|it; K){iu- uMttvUvLvuU. 


s z j ■ 


(64) 


jz IzSz 


The scalar combinations of these matrix elements can be simplified by using the identity 
(UU f ) mim2 = — v^3(lmi; lm2|00) together with the orthogonality relations for Clebsch- 
Gordan coefficients. The resulting formulas are 


{XcJ |'0 t (-|D • cr)xx f (-|D • (t)jJj\xcj) 

~ <5 jo (Xco|'0 t (-|D • <r)xx t (-|D • <x)^|y c0 ), (65a) 

(x c j^\-0ra n xx\-0) m v n MXcj) 
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(65b) 


3 v 2 

{XcM\~0) m o n xx\-0) n a^\ Xc j) 


O (Xco|V ,t (-fD • cr)xx f (-f D • oOV’IXco) , 


(-l) J o (XeolV^-lD • o-)xx f (-|D • cr)^|Xco) • 


(65c) 


Using the relations (|62|) and (|65D , the dominant terms in the decay rates for the XcJ 
reduce to 


r(Xcj) = 


'C& 3 s) 


2 M. 


XcJ 


mt 


{Xco\^crT a x ■ x^Tyixco)^ 


+ 


C'(i- 3 ^)(/i) 


mi 


<Xco|'0 t (—• o’)xX t (-|D • (t)^|x c0 ) 


where (7b-’ 3pj ) is a linear combination of the coefficients C^- ,3p \ n = 1,2,3: 

+ icF p U) + ^c 3 tt* p >. 


In explicit form, these coefficients are 


C (L3po) (/i) = ya^(2m c ) 
C a ’ 3pi) (/i) = ya£(2m c ) 

= y a s(2 m c) 


— In 

27 


0 


4 

27 

— (In 
27 


(i 29 \ a s 

2^ c + ~6~J n/ Y. 

, 4\ cc s 


2m c 3 ' n/ ~ 


— - — (In 
L15 27 


b 


7T 

29 \ a s 


2m c + 15/ 7T 


The constants Qj under the logarithm in these three expressions are 


Qo — 


29 


Q i — 


— 


29 


( 66 ) 


(67) 


(68a) 

(68b) 

(68c) 


(69) 


6 3 15 

Up to corrections of relative order v 2 , the matrix elements in ( J66|) are related to the 
standard matrix elements defined in Ref. |l!] by a simple normalization factor: 


(Xco\^crT a x ■ x'YT'VIxco) ps 4m c (Xco|08( 3 5i)|Xco) , (70) 

(Xco|'0 t (-fD • cr)xx t (-|D ■ <r)^|Xco) « 12m c (Xco\Oi( 3 Po)\Xdo) • (71) 

Using the vacuum-saturation approximation, the matrix element (Xco|^M 3 -fb)|Xco) can be 
expressed in terms of the derivative of the radial wavefunction for the Xco evaluated at the 
origin: 
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(Xco|C > l( 3 -Po)L\co) 




(72) 


27 

2tt 


KcM 2 - 


The corrections are of relative order v 4 . 

We now compare our results for C^- ,3pj ' ) with previous calculations of these coefficients 
to order a^. Our results agree with the Ufa .^ terms obtained in recent calculations by Huang 
and Chao [pj and by Petrelli [jj, who also used dimensional regularization as the infrared 
cutoff for divergences associated with the emission of a soft gluon. Our coefficients need not 
agree with those in the original calculations by Barbieri et al. H||, since they used different 
infrared cutoffs. However any differences in (jA-’ 3p R must be compensated by differences 
in the matrix element (x c0 \'4 > ' <J T a x ■ X 1 crT'VIXeo) in the factorization formula (j66j). This 
implies that the differences between the coefficients C^- 3pjS> obtained with two different 
infrared cutoffs must be independent of J. 

In Ref. |]J, Barbieri et al. calculated the coefficients of the color-singlet matrix element 
in the decay rates of the Xco and Xc 2 using the binding energy of the cc pair to regularize the 
infrared divergences that arise from the emission of a soft gluon. Taking the annihilating cc 
pair to have invariant mass M, the infrared divergence shows up as a logarithm of M 2 — 4m 2 . 
This cutoff can be translated into an equivalent cutoff A on the momentum of the soft gluon 
by replacing ln(4m 2 /(M 2 — 4m 2 )) —> — ln(A/m c ). In Ref. Barbieri et al. calculated the 
decay rates of y c0 , Xci, and y c2 into two-jet configurations defined by an angular resolution 
6 and a fractional energy resolution e. Their result for Xci is independent of <5 and depends 
logarithmically on e. This cutoff can be translated into an equivalent infrared cutoff A on 
the momentum of the soft gluon by replacing ln(2e) —> ln(A/m c ). Expressing the results 
of Barbieri et al. in the same form as in (|68|), the expressions inside the parentheses are 
ln(A/m c ) + Q'j, where the numbers Q'j are 

Q'o = 4 , Qi = i , Qi = h . (73) 


Comparing the coefficients in (|69|) and (|73|) , we see that Qj — Q'j = independent of J. 


Thus the differences between the two calculations correspond simply to different choices for 
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the definition of the color-octet matrix element. The relation between the matrix elements 


in the two calculations is 


(Xco\^<rT a x ■ x^T a i/;\xco) {fl) 


(Xco\^crT a x ■ X^Ty\ X co) {A) 


dim.reg. 


cutoff 



^ (Xco|V' t (-|D • cr)xx t (-fD ■ cry\xco) • (74) 


(1 ^ p\ 

This identification can be verified by repeating the calculation of the coefficient C 2 


in Section 111 using an infrared momentum cutoff A. On the QCD side of the matching 
condition, one must put a cutoff x > A /m c on the integral in (|26|) for n — 0. On the 
NRQCD side, one must put a cutoff k > A on the integral in ([47]). The relation (|74|) is 


identical to the relation between the corresponding production matrix elements defined by 


dimensional regularization and by a momentum cutoff, which was found in Ref. H. 


In Ref. 0 and H. Mangano and Petrclli used the next-to-leading order results of 


Barbieri et al. for XcO and Xc 2 , but they performed an independent calculation for Xci ■ 


Their result corresponds to Q\ = — Since this result disagrees with four other independent 
calculations, we conclude that it is incorrect. 


V. CONCLUSIONS 


In the NRQCD factorization framework, the short-distance coefficients in annihilation 
decay rates and in inclusive production cross sections can be calculated systematically as a 
power series in a s (rn c ). For most applications, next-to-leading order calculations are essen¬ 
tial for accurate predictions. Dimensional regularization is the most convenient method for 
regularizing the infrared and ultraviolet divergences that arise in calculations beyond leading 
order in a s . The generalization of the threshold expansion method to N dimensions allows 
dimensional regularization to be used consistently in quarkonium calculations. In this paper, 
we used this method to calculate the terms proportional to Ti/ctf in the short-distance coeffi¬ 
cients of color-singlet matrix elements in the annihilation decay rates of P-wave charmonium 
states, thus resolving the discrepancies between previous calculations. Our results agree with 
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the original calculations by Barbieri et al. after allowing for differences in the color-octet 
matrix element due to different infrared cutoffs. Our results also agree with recent calcu¬ 
lations using covariant projection methods in conjunction with dimensional regularization. 
In conclusion, the threshold expansion method combined with dimensional regularization 
provides a general and powerful method for carrying out quarkonium calculations beyond 
leading order in a s . 
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APPENDIX A: NONRELATIVISTIC EXPANSION OF SPINORS 


In this Appendix, we present the formulas for the nonrelativistic approximations of 
spinors that are needed to calculate short-distance coefficients for heavy quarkonium de¬ 
cay rates using the threshold expansion method. All the formulas below hold in N spatial 
dimensions. The representation for gamma matrices that is most convenient for carrying 
out the nonrelativistic expansion of a spinor is the Dirac representation: 


^ 1 0 ^ 


/ ,\ 

, y = 

0 a' 


o 

1 


1 

A. 

o 


(Al) 


In the center-of-momentum (CM) frame of the cc pair, their momenta p and p of the c and 
c can be written 


V = {E q , q) , (A2a) 

P = {Eqi -q) > (A2b) 

where E q = Jrn? c + q 2 . The spinors for the c and the c in the CM frame are 
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( 


u(p) = 


v(p) = 


E q + m c 

1 

E 0 + m c 


(Eg + m c ) £ 

q- ° £ 


\ 


\ 

( 

—q • <7 rj 

K (Eg + m c ) r) J 


(A3a) 


(A3b) 


Color and spin quantum numbers on the Dirac spinors and on the Pauli spinors 7] and £ are 
suppressed. If the Pauli spinors are normalized so that rj'ri = = 1, the spinors in ( |A3| ) 

satisfy uu = —vv = 2m c . 

The independent quantities that can be formed by sandwiching 3 or fewer Dirac matrices 
between v(p) and u(p) are 


v(p)u(p) = -2?/ f (q-cr)£, 
v(p)^u(p) = — g™ \2E q rfa^ 


Eg + 7Tl c 


q J rf(q • cr)£ I , 


v(p)(Yl v - YY)u(p) = 4 (^V J - g^gW) I m c 77 + 




Eg + m c 


q> ?7 t (q • cr)^ 


+ g^ {[a>, a% q-cr}Z, 


(A4a) 

(A4b) 

I 

(A4c) 


v(p)(^Yi^ — n t tlr i Vr i ,J ) u (p) 


= „t{[<T>V l R - iT 1 - 


Eg + m c 


Eg + m c 




Eg + m c 


g^{[a\a J ],q- (t}^ 


+ 


4 (<fVY y + g M g ,Jl g 113 + g l,0 g' M g V3 ) (V/h/a^ - q j a . (A4d) 


These expressions can be obtained from the corresponding spinor factors in Appendix A 
of Ref. H by setting LA = —g M = and P^ = 2E q g^° and by taking the hcrmitian 
conjugate. Using the expressions for the spinor factors in (|A4|) , it is easy to carry out their 
nonrelativistic expansions in powers of q. 
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FIGURES 


Fig. 1. The lowest-order cut Feynman diagram in QCD that are associated with the 
process cc —> qq. 

Fig. 2. The lowest-order cut Feynman diagrams in QCD that are associated with the 
process cc —> qqg. 

Fig. 3. The lowest-order cut Feynman diagrams in QCD that are associated with the 
process cc —> gg. 

Fig. 4. The lowest-order Feynman diagram for NRQCD matrix elements such as 
(cc'\'i/j Jf crT a x ■ x^<rT a ^|cc). 

Fig. 5. The next-to-leading-order Feynman diagrams for the NRQCD matrix element 
(cc , |'0^crT a x • x' <J T a il)\cc) when the initial and final cc pairs are in color-singlet states. 
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